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Abstract 

The sine-Gordon model on the half-Hnc with a dynamical boundary introduced by Delius and one of the 
authors is considered at quantum level. Classical boundary conditions associated with classical integrability are 
shown to be preserved at quantum level too. Non-local conserved charges are constructed explicitly in terms of the 
field and boundary operators. We solve the intertwining equation associated with a certain coideal subalgebra of 
Uq{sl2) generated by these non-local charges. The corresponding solution is shown to satisfy quantum boundary 
Yang-Baxter equations. Up to an exact relation between the quantization length of the boundary quantum 
mechanical system and the sine-Gordon coupling constant, we conjecture the soliton/antisoliton reflection matrix 
and bound states reflection matrices. The structure of the boundary state is then considered, and shown to be 
divided in two sectors. Also, depending on the sine-Gordon coupling constant a finite set of boundary bound 
; ^ ' states are identified. Taking the analytic continuation of the coupling, the corresponding boundary sinh-Gordon 
5^ model is briefly discussed. In particular, the particle reflection factor enjoys weak-strong coupling duality. 
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Keywords: Massive boundary integrable field theory; sine-Gordon model; boundary degrees of freedom; refiection 
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1 Introduction 

Two-dimensional quantum field theories with boundary have attracted attention for many years, as they play an 
important role in the analysis of low dimensional statistical systems near criticality or quantum gravity (open 
string approach) . Solving integrable theories restricted on the half-line for various types of boundary is then of 
great interest. Among the known examples, the sine-Gordon model with a non-dynamical boundary has been 
shown to be integrable at classical level [1, 2] as well as at quantum level [3, 4]. This has provided one of the 
simplest examples for which exact results such as factorized scattering theory, boundary spectrum, etc. ..[5, 6, 7, 
8, 9, 10, 11, 12, 13, 14, 15, 16, 17] have been derived. 

^ pascal@yukawa.kyoto-u.ac.jp 
^kkoizumi@cc. kyoto-su.ac.jp 
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Recently [18], an integrable Haniiltonian which describes a sine-Gordon model on the half-line coupled to 
a non-linear oscillator at the boundary has been shown to be integrable at classical level. Using a solution of 
the classical reflection equations (also called classical boundary Yang-Baxter equations) and following Sklyanin 
formalism we have obtained an infinite number of mutually commuting classical conserved quantities, provided 
some specific classical boundary conditions. The existence of such integrals of motion is a sufficient condition which 
ensures classical integrability of the model. Quantum integrability had remained an open question. Although a 
solution of the quantum reflection equations - quantum boundary Yang-Baxter equations (qBYBE) - had been 
proposed for a certain representation of the quantum i?-matrix [19, 20], there was no reason to believe that this 
solution could describe the soliton/antisoliton scattering process bouncing off the boundary. 

Independently, the method of construction of non-local conscirvcd charges in integrable models in the bulk [21] 
has been extended to non-dynamical boundary ones more recently. For instance, the sine-Gordon model [22] and 
the aI^'' afRne Toda theory with imaginary coupling [23] restricted to the half-line have been studied. The existence 
of such non-local conserved charges provides a useful tool to determine the scattering properties of the theory. 
Indeed, the explicit form of the S'-matrix (soliton/antisoliton scattering) and the ii'-matrix (soliton/antisoliton 
scattering on the boundary) is encoded in the minimal solution of quantum Yang-Baxter equations (qYBE)s and 
qBYBEs, respectively. 

In this paper, we study at quantum level the model introduced in [18]. In section 2 we will show that quantum 
boundary conditions take a form similar to the classical boundary ones proposed in [18]. Also, we will construct 
non-local conserved charges corresponding to our dynamical case, which are natural extensions of the known 
ones (non-dynamical) [22, 23] . In our case the non-local conserved charges generate a certain coideal subalgebra 
of Uq{sl2) mixed with the Hciscnbcrg one, where q is the deformation parameter. Further, by specifying the 
representation for Uq{sl2) we solve the corresponding intertwining equations. It provides a new solution Kq{9) to 
the qBYBEs, different from the one proposed in [19, 20, 18]. 

In section 4, we use the minimal solution Kq{6) of these equations'^ to construct the soliton/antisoliton 
reflection matrix KgQ{6) which describes the scattering of the sine-Gordon soliton/antisoliton off the dynamical 
boundary. In particular, taking the commutative limit of the Heisenberg algebra, we check that the ''minimal" 
part reduces to the one obtained by Ghoshal-Zamolodchikov [3] and DeVega-Gonzales-Ruiz [24]. In the generic 
noncommutative case, we show that K^^{9) is a solution of the qBYBEs associated with the sine-Gordon 5- 
matrix. This shows consistency of the nonperturbative analysis based on non-local conserved charges approach. 
To determine uniquely KgQ{0) we impose the boundary unitarity condition as well as the boundary cross-unitarity 
condition proposed in [3] . Studying the singularities of the proposed reflection matrix we identify a finite set of 
boundary bound states. Then, we apply the bootstrap equation to construct explicitly the reflection matrices 
associated with bound states (breathers) scattering off the boundary. 

As we are going to see, with the help of certain algebraic relations between Heisenberg operators the boundary 
unitarity condition is satisfied without specifying the boundary state structure. However, the boundary cross- 
unitarity condition leads to restrict its form in order to preserve certain "physical" principles. We identify the 
corresponding constraints: choosing a continuous representation for the boundary operators difference equations 
are obtained. They restrict the boundary state in two sectors, associated with "even" or "odd" breathers i?„. 

Concluding remarks follow in the last section. There, following Corrigan [12] we use the so-called breather- 
particle identification to obtain the particle refiection factor in the boundary sinh-Gordon model with quantum 
degrees of freedom at the boundary. The result is found to be self-dual under the weak-strong coupling duality 
transformation 6 <-> 2/6. 

2 Classical sine-Gordon field theory with a dynamical boundary 

Let us first recapitulate some main results of [18] without details. There, we considered a Hamiltonian describing 
the interaction between a sine-Gordon field theory restricted to the half-line and a non-linear oscillator living at 
the boundary. The sine-Gordon part describes a relativistic 1-1-1 dimensional self-interacting bosonic field 4>{x,t) 
^These relations, in the dynamical case, are no longer linear in terms of the whole set of operators. 
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with mass m and its conjugate field tt{x, t). The Hamiltonian reads 

HsG = f_Jx(^-n{xf + \{d^<l>{x)f-'^cos(icl>{x)^ (1) 

(cos(/3p/v/2Mo^) e-'^^^W/^ + cos{P ^M^q / 2^2) e'/^^W/^) . 

Expanding the oscillator part for small q and p, these degrees of freedom can be interpreted as position and 
momentum variables, respectively. In this limit, it can be shown that the effective mass of the oscillator depends 
on the value of the field at the boundary and the parameter Mq. 
At fixed time t we introduced the Poisson bracket 

50i 5O2 50i 5O2 



{0^,02} = f 

J — CK 



dx 



dOid02 _ d0^d02 
dp dq dq dp 



-STr{x) S(j)(x) 6(j){x) 6tt(x) 

for any observable Oj . At constant time slices the non- vanishing Poisson brackets in the sine-Gordon field theory 
and for the boundary variables p, q are, respectively 

{7r{x),<P{y)} = 6{x-y) , {p,q} = l (3) 

from which it is straightforward to calculate the equations of motion. However, in order to have the continuity 
of these equations on the left half-line ] — oo,0], the equations for the field 7r{x) leads to the classical boundary 
condition 

9,0(0) = -^^(cosp e-'^'^W/^ - cosg e'^^^/^) . (4) 



where we have introduced the reduced parameters p = "yj^^P and q — ^^^'" g. 

The model with Hamiltonian (1) is integrable: it is possible to construct an infinite number of higher spin 
integrals of motion that are in involution with each other. For instance, assuming the boundary condition above 
(4) the first non-trivial integral of motion reads 

/3 = £ {A^r + 6n{xf{d,^{x)f + {d,^{x))') - ^ {{dM^)f + {dl^{^)f) 

[nixf + 5(5,</.(x))2) cos/3<^(a;) + ^ (cos2/30(x))) dx + 4°™'^-'-^ (5) 



with 



jboundary ^ g3i/3^(0)/2 (^1 coS g + ^ COS^ q] - e'^^(°)/2 f ^ OS p - ^ COS^ qC0Sp+ :^n{Of COS q 

\2 6 / \2 2 2m'' 

^ g-3»w(o)/2 cosp + i cos3 p\ - e-i0mi-^ ^osg - i cos^pcosg + ^<^f cosp 



.2(3 ... . 
— t — smgsmp. 

m 

In general, this Hamiltonian is not real. Such kind of situation is typical in higher rank affine Toda theories 
with imaginary coupling for which the Hamiltonian is also non-hermitian. In these theories, the classical soliton 
solutions are found to be complex. Nevertheless, the energy of these configurations is real [25]. 



3 Quantum sine- Gordon field theory with a quantum boundary 

Let us suppose that there exists a well-defined action describing the sine-Gordon (SG) field theory restricted on 
the half-line coupled with a quantum mechanical system at the boundary. To characterize the dynamics of the 
boundary we introduce, let say, two boundary operators £±{y). We propose the following action 

dy j dxi^—{d^(j)f -2flC0s{pcl))j J dy ^^^^t{y) + Aboundary (6) 
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where the interaction between the field and the diincnsionless boundary quantum operators reads 

and Aboundary IS the kinetic part associated with the boundary operators. Here we used the notation (3 = 13/^/^. 
In the deep UV, this model can be considered as a relevant perturbation of a conformal field theory (for instance 
the free Gaussian field) on the semi-infinite plane with certain boundary conditions at a; = 0. In the limit /x = 
Neumann boundary conditions dx(j)\x=o = preserve integrability. If one turns on the boundary perturbation 
(/z ^ 0), following the arguments of [3] we expect that a quantum analogue of the classical integral of motion 
(5) can be constructed explicitly. Obviously, the boundary operators S±{y) will have to satisfy certain algebraic 
relations in order to preserve integrability at quantum level too. Even if we do not construct explicitly such local 
quantum conserved charges here, analysis of following sections will support integrability of the model. 

There are two alternative description^ of the model (6). On one hand, one can choose the y-direction to be 
the "time" in which case the Hamiltonian Hb contains the boundary contribution, and the Hilbert space Hb 
identified with the half-line y = Const.. Then, correlation functions are calculated over the ground state of Hb, 
denoted |0)b below. 

On the other hand, one can choose the infinite line x = Const, as the "equal time section". The Hamiltonian 
in this case is the same as in the bulk theory on the full line with Hilbert space Ti.. The boundary at a; = 
plays the role of initial condition, and all its information is encoded in the boundary state \B) G H. Correlation 
functions arc written as (0|...|-B) where |0) S W is a (degenerate) ground state of H. Notice that the existence of 
local integrals of motion I^*) , I^**^ for certain values of s implies 

(jW _ =0 . (7) 

In the perturbed boundary conformal field theory (BCFT) approach, the free bosonic fundamental field re- 
stricted to the half- line can be written in terms of its holomorphic/anti-holomorphic components 

(l){x,y) = ip{z) + Tp{z) . 

If we denote the expectation value in the BCFT with Neumann boundary conditions (...)oj then the holomorphic 
components are normalized such that 

{ip{z)ip{w))o = -21n(z - w), {if{z)(f{w))o = -21n(z - w), {ip{z)(p{w))o = -21n(2 - w) . 

To find the boundary conditions at quantum level, we can consider the expectation value of the local field 9j;0(O, y) 
with any other local field in first order of boundary conformal perturbation theory ^ 0. Due to the presence of 
the boundary operators £±{y), we write the ground state |0)b as 

|0)b = \vac)B ® \Q)bcft + 0{ii) , (8) 

where \Q)bcft € Hbcft and \vac)B denotes the effective boundary ground state which satisfies the 
Schrodinger equation 

H{£+, £-)\vac)B = E^^i\vac)B (9) 

with effective ground state ^ energy E^-^^ . Here, H {£+,£-) is an effective boundary Hamiltonian. In this ap- 
proximation, following the analysis of [23] with eqs. (8) it is straightforward to show that, in first order of 
boundary conformal perturbation theory, the quantum boundary condition takes a form similar to the classical 

^As we will sec later, integrability of the QFT (6) requires that one can not turn on the boundary perturbation independently of 

/i. 

^It seems important to us to recall the analysis of [3]. 

^Taking the perturbative limit /3 — > and using an explicit realization of the boundary operators in terms of Heisenberg operators 
(see further sections), one can show that the "pure" boundary wave function satisfies certain difference equations. 
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one. Furthermore, using scaling arguments ^ it is possible to show that this condition is satisfied at all order in 
perturbation theory. It reads 

B(O|9,0(O,y)...|O)B -^27r/3M^/2 B(0|(£-(y)e*'^^(°''')/^ - f+(y)e-''^^(°'^)/2)...|0)B • (10) 

The method of constructing non-local conserved charges in the model (6) follows the line presented in [23]. So, 
we refer the reader to this paper for more details. Here, the main difference is that S±{y) are operators. If we 
denote the "bulk" non-local charges [21] 

Q± = ^p dx{J± -H±), '5± = ^ dx{J± - H±) , 

they can be expressed in terms of holomorphic/anti-holomorphic part of vertex operators as^ 

J± =: exp(±-^i^) : and H± = —Awfj,^ — - : exp |^±i — /3 ) =F ifDip 
J± =: exp(=F-^(^) : and H± = —A-Kji-^ — ^ : exp y-^ — j3j (p^ i(3(p 

Together with the "bulk" topological charge T^uik = 4k I^oo ^^^Y generate the quantum enveloping 

algebra Uq{sl2). In the theory on the half-line, those charges are no longer conserved. Instead, using (8) and the 
method of [22, 23] it is straightforward to show that the two combinations 

Q± = Q± + Qt + £±{y)Q'^^ with ^ = ^/ '^^^^'^ 



CSO 



and ^2 

£±{y)=n'/'-\^£±{y) (12) 
1 — 



are conserved to all orders in boundary conformal perturbation theory framework. Here, the deformation param- 
eter is defined by 

q = exp(-2i7r//3^) . (13) 

As the sine-Gordon model possesses a single two dimensional soliton multiplet denoted by |V'±(^))) we are now 
interested in two dimensional representations of Uq{sl2) on asymptotic soliton states, denoted irg below. Here, 6 
denotes the rapidity of the soliton/ antisoliton. To determine the meaner in which non-local conserved charges (11) 
are represented on asymptotic single-soliton states, we use the notation of [21]. Furthermore, in the asymptotic 
limit y — > ±oo on the boundary we will assume 

S+{y = ±oo) = S+ and B-iy = ±oo) = £- . (14) 

In the first Hamiltonian picture y is identified with the time axis. In other words, the condition (14) means that 
boundary operators lead to the same operators in the far past and future. We must keep in mind that are 
operators satisfying [£±,g\ = Q and 7T${£±g) = £±n0{g), Vgi e f7q(s/2)- Thus we have 

MQ±)+ = £±1^' , MQ±)- = ce±^' , MQ±)+ = ce^^' , MQ±)- = £±q^' 

with 

\ = 2/i3^-l and ^ i2^i{q'^ - 1)/X^ . (15) 



^Dimensions of both side must be equal. Due to the form of the perturbing operator one finds that the only term that can appear 
on the right hand side must be linear in the parameter fi. All other terms must vanishes. 

*On the contrary to [23], here we denote the fundamental field one the whole line (l>{x,y) and the one restricted to the half-line 
(p{x, y). Then, the chiral components are related in the following way: <p{x, y) = '}>{x, y) + lp{—x, y), ^{x, y) = Tp{x, y) + <fi{—x, y). 



5 



For the remaining symmetry algebra of the QFT (6) generated by the non-local conserved charges given above 
(11) we are looking for a solution to the intertwining equation 

Koi{0)MQ±)l = ^-e(Q±)!i^oc(^) (16) 

where indices {S,u,Q G {±} refer to two dimensional representations of Uq{sl2) on asymptotic soliton states 
\tp±{0))- Then, using this representation the solution Ko{9) is written as a 2 x 2 matrix with entries expressed 
in terms of the boundary operators. In the non-dynamical case, the entries are just analytic functions of 6 as £± 
are c— numbers. For further convenience, let us define 

Ko+{e) = A{e) , Kot{e) = b{9) , 

Ko+ie) = D{e), KoZie) = Eie). (17) 

After some calculations, we find that the entries of the minimal solution Kq(6) of the intertwining property (16) 
takes the following form: 

A{e) = {q-'e^'£+ - qe-^'£^){q - q'^)/2c , 
E{9) = {q-'e^'£_ - qe-^'e+){q - q-')/2c , 

B{e) = ( - c^g-ig^^^ - c\e-^^' + i^ili(g-if_£+ - qS+e-))/2c^ , 

D{e) = ( - c^-'e"^' - c\e-^^' + ^-^{-q£_£+ + q-'£+£_)) /2c^ , (18) 



if the boundary operators satisfy certain commutation relations with the entries B{6) and D(6). Among these, 
we obtain the commutation relation ^ 

[£+£.,£.£+] = c\q + q-'f{£l - f!) . (19) 

Here, we propose a realization of the boundary operators at y = ±oo in terms of operators {p, q\ which belong 
to the Heisenberg algebra [j5, g] ~ a in order to have (19) satisfied. Due to (14) we denote p(±oo) = p and 
g(±(x)) = q. Indeed, one can check that setting 

£4. = ±2ei7y coshp and = ±2er/\/ coshfj where ±<7 = cxp(q/2) (20) 

with the normalization e^y = —c'^/{q — q~^)^, all the equations (16) can be satisfied for a certain relation between 
the SG coupling constant and the boundary quantization length a. However, according to eqs. (20) there is an 
ambiguity in the definition of the sign of the boundary perturbation. Depending on each sign, the boundary 
quantization length is fixed to 

(g'2 _ I) (02 _ 2) 
[p,q] = a mod {Am) with a = iA-K^^ — - (+) or a = i2-K-^^ (-) . (21) 

In the model (6) with (12) and (20), euv is fixed by the "bulk" non-local conserved charges algebraic structure 
which gives 

2 i2{l-q^) 

To relate the UV parameter euv and the IR soliton mass parameter M we can use their exact relation calculated 
in [26]: 



r(/3V2) r M^T {1/2 + 1/2A) 



^ ~ ^r((2 - /32 )/2) L 2r(i/2A) J • ^ 

9For generic values of the coupling and explicit boundary operators dependence in the boundary scattering description, this 
condition is necessary. Also, we do not assume specific values of the coupling here (e.g. free fermion point, reflectionless point), for 
which new solutions may be constructed. 

^"Notice that this solution also works, obviously, if one changes a — * —a. 
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It follows that the strength of the boundary perturbation is fixed by the bulk soliton mass M . This is consistent 
with the classical model (1) in the sense that the only free parameter is the mass m^^. In other words, choosing an 
arbitrary boundary mass would destroy intcgrability. Physically, this phenomena is not so surprising: to preserve 
integrability the effects brought by the SG model and the non-linear oscillator (energy transfer, for instance) have 
to compensate each other, which is characterized by the exact relation between the quantization length a and the 
SG coupling constant in (21). 

At specific values of the SG coupling constant 2/(n + 1) with n e N the operators £± commute, and we 
have = 1. Then, there is no need to assume some special relation involving c and the boundary parameters like 
(19), and £± remain free. Consequently, Ko{6) is reduced to a 2 x 2 matrix where £± are c— numbers. Expanding 
Ko{6) in powers of {q^ — 1) we obtain 

Koi9) ~ 2cosh(2A^)P/(g - q-') + iKZ^-dyJO) , 

where P is the 2x2 permutation matrix. The first (singular) term is a trivial solution of the intertwining 
property (16) whereas the second term is the minimal reflection matrix obtained by Ghoshal-Zamolodchikov and 
DeVega-Gonzales-Ruiz [3, 24]. 



4 Boundary Yang-Baxter equations and factorized scattering theory 

For the coupling constant < 2, the bulk sine-Gordon model in 1-1-1 Minkowski space-time is massive and 
intcgrable. The particle spectrum consists of a soliton/antisoliton pair ('0+(^?),'0_(0)) with mass M and neutral 
particles, called "breathers", Bn{9) n = 1, 2, < A. As usual, E = M cosh9 and P = Msinh6' the energy and 
momentum, respectively, of the soliton/antisoliton. In this model, H is the Fock space of multiparticle states. A 
general A''-particles state is generated by the "particle creation operators" A^. 

\AaAOl)...Aa^{ON)) = AaAOl)...Aa^{ONm (24) 

where characterizes the type of particle. The commutation relations between these operators are determined 
by the ^-matrix elements, which describe the factorized scattering theory [27]. Integrability imposes strong 
constraints on the system which implies that the 5-matrix has to satisfy the quantum Yang-Baxter equations. 
For instance, the SG soliton/antisoliton scattering 5- matrix can be written as [27] 

SsG{e) = R{6)p{-ie)/i , (25) 

where we introduced the four dimensional representation of the trigonometric solution of the quantum Yang- 
Baxter equation written in terms of the standard Pauli matrices CTfe, A; = 1, 2, 3: 

R{9) = ^ (/ (8) / -h 0-3 era) + ^(10/- 0-3 00-3) + ^(cti (8)CTi +£72 «)fT2) (26) 

with 



a{e) = sinh(a7r - XO) , b{e) = sinh(A6l) , c{e) = sinh(iA7r) . 

The amplitudes b{6) and c{9) possess simple poles in the physical strip < 9 < in located at 9n = in — inn/X for 
n = 1, 2, ... < A. They are interpreted as the neutral bound state (breather) B„. The factor p{u) in (25) ensures 
unitarity and crossing symmetry of the ^-matrix. It reads 

piu) = --nxni - xuMni - x + xum n ^'^^'Lvr^ ' (^7) 

TT ^J-^ Fi{Q)Fi{'k) 

^^Whatever the mass of the eflfective oscillator is, the relevant quantity is its frequency a; = m/2 [18]. Using the particle-breather 
identification, one can relate m associated with the lightest breather with the mass m of the fundamental particle. 
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with 

r(2/A - Au/7r)r(l + 21X - Xu/tt) 



Fi{u) 



T{{21 + 1)A - Xu/Tr)T{l + {21 - 1)A - Xu/tt) 



Using the bootstrap equation, the ampHtudes associated with (anti)soUton-breather tjj±Bn and breather-breather 
BnBm scattering have been calculated in [27]. 



4.1 Soliton reflection matrix for the dynamical boundary 

Taking the first Hamiltonian picture, we can consider y — it in the QFT (6). If one uses the approach of [3], 
asymptotic states are now generated using Aa{9) acting on the ground state |0)s GHb- For boundary integrable 
field theories, the action of the creation operator on the boundary is characterized by the reflection matrix. It 
describes the scattering of the soliton/antisoliton on the boundary and is constrained by the so-called reflection 
equations, i.e. qBYBEs for a certain choice of representation. Depending of the quantization condition in (21) 
which determines the sign e e {±} of the boundary operators (20) - the sine-Gordon reflection matrix KgQ{6) 
must satisfy the qBYBEs 

1 2 

where we used the notations K = K^I, K = I ® K. Recalling the definition of the entries (17) we have fourteen 
functional equations: 



ii) 


a-c+ {BD' - 


- B'D)+a-a+[A,A']=0 , 








(u) 


b-b+[A,E'] 


+ c^c+[E,E']+ c^a+{DB' 


- D'B) = , 






(Hi) 


c_b+{EA' - 


E'A) +b_c+{AA' - E'E) H 


- b_a+[B,D' 


1=0, 




(iv) 


b^b+AD' + , 


C-C+ED' + c^a+DA' - a_ 


a+D'A- a_ 


c+E'D 


= 0, 


(v) 


b_b+B'A + 


c_c+B'E+ c_a+A'B- a_ 


Ai+AB' - a_ 


c+BE' 


= , 


(vi) 


b-b+D'E + 


c-c+D'A+ c-a+E'D- a. 


-a+ED' - a. 


.c+DA' 


= 0, 


(vii) 


b-b+EB' + 


c-c+AB' + c-a+BE' - a. 


-a+B'E- a- 


c+A'B 


= 0, 


{via) 


b-a+BE' + 


c-b+EB' + b-c+AB' - a_ 


-b+E'B = , 






{ix) 


b-a+A'B + 


C-b+B'A+ b-C+B'E- a_ 


.b+BA' = , 






(x) 


b-a+E'D + 


c-b+D'E+ b-c+D'A- a. 


-b+DE' = , 






{xi) 


b-a+DA' + 


C-b+AD' + b-C+ED' - a. 


.b+A'D = , 







where we used the shorthand notations a_ = a{6 — 6'), a+ = a{9 + 6') and similarly for b and c as well as A = A{6) 

and A' = A{9') and similarly for B,D and E. The remaining three equations arc obtained from (i), (m), (Hi) 
through the substitutions A E and B ^ D. Straightforward calculations show that K^q{9) takes a form 
similar to Ko{9). Indeed, the reflection matrix in the dynamical boundary case can be written 

Ki%{9)=Kl,'\9)Yi9), (29) 

where the function Y{9) has to be determined using some physical assumptions (sec below). In the two dimensional 
representation of Uq{sl2), the matrix Kq"^"^ is obtained from (18) with the following substitutions^^: 

f+ ecoshp , f_ ecoshg , q — > e"^"^ and c — > sin(7rA) . (30) 

Here, the quantization condition is given in (21). In [3], Ghoshal and Zamolodchikov proposed the use of the 
"boundary unitarity" and "boundary cross- unitarity" conditions to determine the overall factor Y{9) associated 

^■^Notice that for specific values of the SG couphng constant such that g is a root of unity, it is easy to construct a finite dimensional 
representation for the operators A, B, D, E in the reflection matrix. 
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with a non-dynamical boundary. Then, let us first consider the boundary unitarity condition. Independently of 
any representation for the boundary ground state |0)b, in our case it is natural to consider the combination of 
operators 

4CW4'^6(-^) = ^6I with {a,b,c}G{±} (31) 

as a generalization to the dynamical boundary case of the condition associated with the non-dynamical one. Here, 
the operator I denotes the identity which acts trivially on the boundary ground state. In terms of the entries 
defined in (17) these relations becomes 

A{9)A{-9) + B{e)D{-9) = I/Y{e)Y{-9) , 
D{e)B{-9) + E{9)E{-9) = I/Y{9)Y{-9) , 
A{9)B{-e) + B{9)E{-9) = , 
D{9)A{-e) + E{e)D{-9) = . 

As the entries are noncommuting betwcx^n each others, in general one would expect the product (31) to be ill- 
dcfincd. However, due to the specific form of the boundary operators one can show that changing 9 —9 does 
not affect the relations above. Also, the first two relations give the same result. Then, using algebraic relations 
like 

(g2 + q-^)£+£_£+ - e\£_ - e_el = -c'iq + q-^f£- , 
{q^ + q-'^)£-£+£- - £+£l - £H+ = -C^{q + q-^f£+ 

after some calculations we arrive at the condition 

Y{9)Y{-e) = [sinh2(2A6l) - smh^{inX)] . (32) 

Compared to the non-dynamical boundary sine-Gordon, the situation here is rather simple: there are no free 
parameters left in the unitarity condition (32) appart from the coupling constant which appears through A. 

Let us now consider the boundary cross-unitarity condition. In the alternative Hamiltonian description, x 
is now interpreted as "time" and the space of states is the same as in the bulk theory, i.e. H. The "initial" 
(boundary) condition at a; = is encoded in \B). If the theory is integrable (7) must be satisfied. It follows that 
|B) is a superposition of asymptotic states of the bulk theory constituted by pairs of particles of equal mass but 
opposite rapidities [3]. Thus, we define the boundary state as 

|B)=AA(|Bo) + ^<?"i340)|B„) + i r^deK'^'^l{z^-9)Aa{9)A,{-9)\B,,) + ..) . (33) 

Here A/' is a normalization coefficient and the states |Bai...aiv) where Vg is a representation for bound- 

ary operators. The second term corresponds to the contribution from the zero-momentum particles, namely 
the breathers. The coefficients indicate their contribution. Considering in particular the soliton scattering, 
arguments of [3] can be applied and we assume the boundary cross-unitarity condition 

K^sllii^f^ - ^) = SsG% (2^)4'^^-I' {in/2 + 9) , (34) 

as the "in" and "owt" states are related through the 5-matrix. Notice that this relation is linear in the boundary 
operators. For simplicity, let us introduce two meromorphic functions Yq{9) and Y\{9) which solve the functional 
equations (w = —i9): 

Yo{9)Yo{-9) = 1 and YQ{m/2 -9)= sin(A(7r - 2u))p{2u)Yo{m /2 + 9) ; 

Yi{9)Yi{-9) = [sinh^(2A6l) - smh^{inX)] and Fi(z7r/2 -9) = Yi{in/2 + 9) . 
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Then, it is not difficult to show that for the clioice 

Y{e)=Ya{e)Y^{9) 



(35) 



in (29) , the boundary unitarity and boundary cross-unitarity conditions are satisfied. Following Ghoshal-Zamolodchikov 
notations we finally obtain 



where we used [3] 

oo 

= n [ 

and we introduced 



Yo{0) = Ro{u)Go{u) 



T{AXk - 2Au/7r)r(l + 4A(fc - 1) - 2Xu/n) 
L Lr(A(4fc - 3) - 2Au/7r)r(l + A(4fc - 1) - 2Xu/'k) 



(36) 



/{u -u) 



Go{u) = n [ 



fe=i 



r(l + (4fc - 2)A - 2Aw/7r)r((4A: - 2)A - 2AM/7r) 
r((4fc - 4)A - 2Au/7r)r(l + 4fcA - 2Xu/Tr) 



/{u -u) 



Notice that the factor Ro{u) contains poles in the "physical strip" < m < 7r/2 located at u„ = n7r/2A for 
n = 1, ... < A . They are associated with zero-momentum soliton-antisoliton bound states (see last section). Also, 
we find 



cr(7r/2 + 7rA/2, u)cr(7rA/2, u) 
sin(7rA) 



(37) 



where [3] 
(7{x,v) -- 



cos X 



^ rr(l/2 + {21 - 1)A + x/tt - A?./7r)r(l/2 + {21 - 1)A - x/n - Xv/n) 



cos{x + Xv) fj^ L r(l/2 + {21 - 2)A - x/n - Xv/tt)T{1/2 + 21X + x/n - Xv/ir) 



-V) 



satisfies the relations 

a{x, v)a{x, —v) = cos^ x/ (cos(a; + Xv) cos{x — Xv)) , 



a{x, 7r/2 — v) = a{x, 7r/2 + v) 



Prom (29) it is clear that there are no resonance states in comparison with the non-dynamical SG model. However, 
the factors a{n/2 + 7rA/2, u) and cr(7rA/2, u) brings an infinite number of singularities at real values of u. In the 
"physical strip" < u < 7r/2, depending on the value of the SG coupling constant, they are located at: 



even 



T 



for 



< A; < 



A 



and 



TT {2k - l)7r 
2 2A 



for 



< fc < 



A + 1 



Denoting Eq as the ground state energy, wc identify these poles as boundary bound states with energy E^"^ —Eq = 
M cos{u°^'^ ). For any k, we conclude that the boundary bound states are stable. 



4.2 Bound state reflection matrix for the dynamical boundary 

Above, we considered the reflection matrix of the soliton/antisoliton off the boundary. We now turn to the 

reflection matrix of the boimd states (breathers) Bn{9). To calculate it, we use the boundary bootstrap method 
introduced in [3, 5] and mainly follow [4]. Using "particle creation" operator formalism, the reflection matrix 
defines the commutation relation 

Bn{9)B = b!-;;^ B,,{-e)B 

where Bn{9) creates the bound state i?„ with rapidity 9. Then, the boundary bootstrap equation reads [3] 

/m.4C (0 + 0n/2)SsG%i {29)K^;li {9 - 9J2) = f^^^R^^^ {9) (38) 
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where /j"^^ are vertices satisfying /^_(— 1)" = /"_,_ and f!^j_ = 0. Furthermore, the solution of (38) has to 
satisfy the boundary unitarity and boundary cross-unitarity conditions 

i?^"'(zV2 -6)= i?^'^(z7r/2 + 0)5("'")(20) , (39) 

respectively. Here, 5*^"'"' (26*) denotes the scattering amplitude for B^ + B^ B^ + B^ process [27]. In particular, 
it possesses a pole in the physical strip located at ^ = mr/2X. It is straightforward to show that 

Yiie + 0j2)Yiie-0j2) = ,,,,(2,,) 3inh(2A. - 2^X.) ^^"^ ^^/^ + ^^2, u)Si-^ (.A/2, .) . 
Here we used the notations of [15, 16] i.e. 

<^M = 4fe#n^#^^ and 5(")(.,.)=nlf-i±^| 

where we have the shorthand notation (.x) = sin(u/2 + a;7r/2)/ sin(M/2 — xTr/2). Taking into account the non- 
diagonal part of the soliton reflection matrix, we finally obtain after some calculations the result 

R^^\0) = (-l)"i?(")(M)5(")(0,u)5(")(7r/2,u)5(")(7rA/2 + 7r/2,u)S(")(7rA/2,u) (40) 

which, clearly, does not depend on the boundary operators £± . As expected, the factor i?Q"-' (u) contains the poles 
in the physical strip located at u„ = 7r/2 — mr/2X for A > 1. They are associated with the breather-breather 
bound states i?2n- 



4.3 Restriction of the boundary state 

Up to now, we only considered the scattering properties without imposing any restriction beyond (7) to the 

boundary state \B). However, it is necessary to exclude those which would violate certain physical principles, i.e. 
it must be in accordance with the SG scattering processes. For instance, the term Ro{—i9) in eq. (35) contains 
poles in the "physical strip" < < iTr/2 located at = m7r/2A; n = 1, ... < A. They are associated with 
the zero- momentum breathers i?„(0) that should appear as bound states of soliton- antisoliton. At these values of 
the rapidity any off-diagonal entry of (29) which appears in (33) applied to a state |B„) will not, in general, give 
a null result. As breathers are neutral, the state |B„) must however be a null vector of B{0n) and D{0n), i.e it 
must obey the condition 

B((?„)|B„)=0 and D{0^)\B^)=O (41) 

for all n = 1, 2, ... < A. In the non-dynamical boundary case [3] this was trivially satisfied due to the form of the 
off-diagonal entries. Combining these two conditions we obtain 

[coshp, coshq] |B„) = and coshp cosh(7|B„) = (-1)" cos^(7rA)|B„) . (42) 

To characterize each state |B„) explicitly, we introduce the one dimensional "space" representation {|x)} for the 
Heisenberg operators defined by 

p\x) = -2iTrXd/dx , q\x) = x\x) and (x|B„) iE *„(-jx/7r) . (43) 

From the first equation in (42), it is straightforward to obtain a difference equation satisfied by 4'rs(y) for y = 
—ix/iT. Obviously, at that stage n does not appear explicitly. It roads 

^n(y + 2A) ^ sin(7r(y - A)) 

*„(y-2A) sin(7r(y + A)) ' ^ ' 
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Up to an irrelevant overall eonstant, it can be solved by the integral representation 



*„(y) = exp / 
Jo 



dt sinh^ {{y -l/2)t) 



for 7^e(y) < 1 + A (45) 



t sinh(i) cosli(2At) 

and is defined through analytic continuation otherwise. The second equation will fix the "allowed" values of y. It 
reads 

(-l)"2cos2(^A)«'„(y) =cos(7ry)(4'„(y-2A) + «'„(y + 2A)) for y S {y^")} . (46) 

The condition (46) is divided into two sectors, "odd" n breathers and "even" n breathers. Then, using the 
definition (45) the boundary state |B„) can be decomposed on the finite set of elementary "odd" and "even" 
states n e {2p, 2p+ 1} with |B„) = N„^'„(y("^)|y(")) (N„ are some normalization constants) which, due to (33), 
implies that ii'^^^(6'„)|B„) will be automatically block diagonal. 

Notice that in the commutative limit A = in (44) the function ^n{y) remains arbitrary and (46) leads to 
y = n. In other words, the boundary state \B) is not restricted anymore in agreement with [3]. 



5 Concluding remarks 

As was suggested in [3] , boundary integrable quantum field theories including degrees of freedom at the boundary 
can be constructed. Here, we provide such an example corresponding to a sine-Gordon quantum field theory 
coupled to a non-linear quantum oscillator at the boundary, which follows the work initiated in [18]. For certain 
boundary conditions, integr ability of the model is preserved at quantum level for an exact relation between the 
sine-Gordon coupling constant and the boundary quantization length. In particular, this later quantity fixes 
the overall sign of the boundary perturbation. The corresponding soliton/antisoliton refiection scattering matrix 
-^SG (^) constructed explicitly using an extension to the dynamical boundary case of the method based on 
boundary non-local conserved charges [22, 23]. Bound state reflection matrices are also constructed explicitly. 
Using some physical constraints, the general form of the boundary state is shown to be restricted. 

To study this model further, extensions to the dynamical case of known nonperturbative approaches could be 
useful. In the non-dynamical case, they have provided an efficient way of studying boundary effects in finite-size 
system. For instance, TBA analysis [9] or boundary reflection amplitude method [28] from which the effective 
central charge and boundary ground state energy can be deduced. Comparison of both approaches provides 
important checks of the exact relation between UV and IR parameters as well as the boundary ground state 
energy. Here, such methods extended to the dynamical case would obviously provide important information. 

Due to the form of our model, one may wonder if it could not be related with certain massive extension of the 
fleld theory description of the Kondo model [29]. In [30] two massive versions of the anisotropic spin 1/2 Kondo 
model have been proposed. It is thus interesting to compare results from [30] (see also [31]) to ours. The first 
massive Kondo (MK) model is one natural generalization of the integrable massless Kondo model for arbitrary 
spin studied in [29]. The MK Hamiltonian reads 

Hmk = \f dx{{-K{x)f + {d,cj){x)f - Gcos(/3<^(x))) + AolS+e*'^*^")/^ ^ s_^-^0m/2^^ (47) 

where S± form a spin 1/2 representation of the go— deformed quantum algebra 5'Ugo(2). The sUqQ{2) relation are 

s —S 

[5„ S±] = ±2S± , [S+, S-] = and go = exp(«/3V8) . (48) 

qo - % 

This model has been studied in details in [31] at the free fermion point 0^ = Air and in [30] at the reflectionless 
points for which it was claimed to be integrable. Unlike the massless version, it is not integrable for arbitrary j3. 
On the other hand, in [30] the modified massive Kondo model (mMK) with Hamiltonian 

jjmMK ^ ^J°dx(^{TT{x)f + {d,^{x)f-GcOs{P^{x))) 

-iXo{S+ cos{P{cf>{0) - 4>)/2) - S_ cos(/3(<^(0) - <^*)/2)) (49) 
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was proposed. Here = — ^(f — i4'o) with 0o a free parameter and the couphng Aq is related to the bulk coupling 
G. This model is believed to be integrable for arbitrary /3 and a reflection matrix has been conjectured. 

As was mentioned in [31, 30], it is important to notice that the analysis supposed that the boundary scattering 
does not involve boundary operators. Here wo didn't assume this condition, which is crucial for our analysis and 
leads to a new solution (29) of the boundary Yang-Baxter equations. Then, it would be interesting to see in such 
case if more solutions can be obtained if a "simpler" bidk scattering is considered instead (reflectionless points). 
Notice that the free fermion point and the reflectionless points are special cases for which our operators (20) 
£± commute, i.e. the Kq reflection matrix reduces to Ghoshal-Zamolodchikov and DeVega-Gonzales-Ruiz [3, 24] 
result (see end of section 3). 

Other remarks can be made. First, the model niMK is integrable for arbitrary /? only for a special relation 
between G and Aq, and the deformation parameter go is expressed in terms of the SG coupling constant (3. It 
is also the case for our model. Secondly, in the MK and mMK models the boundary dependence only occurs 
in the ODD factors. With the boundary dependence isolated, the massless limit follows by taking the limit 
9 — > oo, which leads to the (massless) Kondo model reflection matrix. However, in our model the boundary 
dependence occurs through the operators cosh(p), cosh(g) in the minimal (non-CDD part) of the reflection matrix 
(29). Consequently, taking the massless limit of our model requires an analysis which goes beyond the scope of 
this paper. Finally, for q a root of unity it is easy to construct a flnite-dimensional cyclic representation for £±. 
This property might have a certain interest: partition functions of the massless boundary sine-Gordon and the 
anisotropic Kondo model can be related if the boundary spin is in a cyclic representation [29] . 

Among other models, it is now quite natural to consider the analytic continuation of the boundary sine-Gordon 
model studied above. It should provide one of the simplest boundary integrable QFT, the sinh-Gordon model 
restricted to .t < with quantum degrees of freedom at the boundary. Its spectrum consists of only one particle 
with mass m and the bulk scattering ^-matrix for a pair of particles reads [32] 



( - l) 26^ 



In particular, this amplitude is self-dual under the weak-strong coupling duality transformation 6 <-> 2/6. Notice 
that the bulk ^-matrices associated with higher simply laced Toda theories enjoy the same property. As was shown 
in [12, 13, 14] for non-dynamical integrable boundary conditions, the sinh-Gordon reflection factor obtained from 
analytic continuation of Ghoshal's result [4] is not self-dual. 

Here, taking the analytic continuation $ ^ ih in (40) for n = 1, the fundamental particle reflection factor 
takes the simple form 



KshG{e) 



4 J 



Then, up to an overall sign in front of the boundary perturbation, integrability is ensured for the boundary 
quantization length (after changing /3 — > i6 in (21) (— )) fixed to 

[p.,q\^ashG with ashG=i^'^/B. (50) 

Consequently, weak-strong coupling duality of the boundary integrable sinh-Gordon model corresponds to the 
boundary operator-coupling transformation 

{P*A\ b*} = {hp/V2, hq/V2, 2/h} . (51) 

Although we didn't discuss them here, it seems to us that our analysis can provide a starting point concerning 
certain open problems: 

• A different choice of trigonometric i?-matrix gives the Boltzmann weights of the six-vertex model. It is known 
to be related to the XXZ and the XXX (in its rational limit) spin chains. Applying our result to this case, the 
value of the boundary quantization length will define the value of the anisotropy which preserves integrability of 
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the model. Consequently, we hope our solution Kq{9) which satisfies the corresponding qBYBEs can be useful in 
the study of spin chains with dynamical boundary conditions 

• Statistical models with extended line of defects have attracted attention as such inhomogeneities affect the 
critical properties of the pure statistical systems. In the continuum limit, the scattering theory of their massive 
excitations is described in terms of the bulk scattering amplitudes as well as those associated with the interaction 
of particles with the defect line. If integrability is preserved, they reduce to reflection-transmission amplitudes 
which satisfy the so-called reflection-transmission equations [35]. Similarly to [18] and in the present work, one 
may be interested in solutions of these equations with quantum boundary degrees of freedom in the defect, taking 
a reflection matrix of the form (17), (18). We will discuss this problem elsewhere [36]. 
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